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Abstract. We present here an efficient method which systematically reduces the rank 
of the augmented space and thereby helps to implement augmented space recursion 
for any real calculation. Our method is based on the symmetry of the Hamiltonian in 
the augmented space and keeping recursion basis vectors in the irreducible subspace 
of the Hilbert space. 



PACS numbers: 71.20,71.20c 
1. Introduction 

In a series of earlier communications Pl Ej we have proposed the augmented 
space recursion (ASR) as an efficient computational technique for the calculation of 
configuration averaged electronic properties of random binary alloys. The method is a 
combination of the augmented space formalism (3j and the recursion method of Haydock 
et al When coupled with the local density functional approximation based tight- 
binding linearized muffin-tin orbital 1^ it provides a ffist-principles determination 
of electronic structure and total energy calculations for disordered alloys. Effects 
such as short-ranged ordering and local lattice distortions due to size mismatch of 
the constituent atoms can very easily be incorporated in the methodology [3 IHl El- 
Recently the method has been applied to study dispersion and line widths for phonons 
in disordered alloys. Since recursion can take into account disorder in beyond nearest- 
neighbour force constants as well as large environmental effects, this application was a 
step forward in the theory of phonons in disordered alloys [10] • It has been pointed out 
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repeatedly in the above publications that the ASR with its terminator approximation 
[m IT^ goes beyond the standard single-site mean field theories and accurately takes 
off-diagonal disorder into account. A number of attempts have been made for practical 
implementation of mean-field theories beyond the single site approximation. One of 
the most sophisticated developments among the mean-field theories, is the travelling 
cluster approximation of Gray and Kaplan IT^ and its latest version : the itinerant 
coherent potential approximation (ICPA) ^Sj, where a small number of classes of 
selected excitations in augmented space are preserved. In spite of the simplification 
so introduced, the resulting equations are still formidable enough to discourage further 
generalization to larger clusters beyond the nearest neighbour. In the cluster CPA 
proposed by Mookerjee JOl EI, in which one retains disorder fiuctuations only from 
a chosen cluster and replaces the rest of the system by an effective medium, one is 
restricted by the limitation of solving number of coupled self-consistent equations. The 
number of equations increases with the size of chosen cluster, which again discouraged 
extension to larger clusters beyond the nearest neighbour. This brings us back to the 
attempt at the implementation of the ASR for large environment fluctuations, implying 
accurate evaluation of a large sequence of recursion coefficients. We shall discuss this in 
some detail in the following section. 

The ASR, in its practical implementation to real alloy systems, has several practical 
problems. The formalism maps a disordered Hamiltonian, described in a Hilbert space 
spanned by a countable basis in which preferably its representation is sparse, onto 
an ordered Hamiltonian in an enlarged space. This enlarged or augmented space is 
constructed by taking the outer product of the conflguration space of the random 
variables of the disordered Hamiltonian together with original Hilbert space in which 
the disordered Hamiltonian was described. The main difficulty in the implementation 
of recursion in this augmented space is its enormous rank. A system of N sites with 
binary distribution has an augmented space of rank N x2'^ . This enormous rank of the 
augmented space and its description on the computer is what had discouraged earlier 
workers in implementing the ASR, although the augmented space formalism was first 
proposed about thirty years ago [3j. The second problem was the recursive propagation 
of errors. Very small errors during recursion (inevitable in numerical implementations) 
tends to add up leading to loss of orthogonality in the recursion basis. This leads 
to appearance of ghost states. Haydock and Te PH] has overcome this difficulty in 
proposing the dynamical recursion. We refer the reader to their paper for details. 

In this communication we shall describe the reduction of the effective rank of the 
augmented space on which the recursion is carried out, using the local point-group 
symmetries in both the Hilbert space and the configuration space. Such reduction is 
well known in various manipulations in reciprocal space for ordered systems. However, 
such a reduction has not been implemented in real space recursion, although it had 
been proposed by Gallagher jTHI quite some time ago. It had been shown by Gallagher 
that if we start recursion with a state belonging to an irreducible subspace of a Hilbert 



Symmetry reduction in the TB-LMTO-ASR for random binary alloys 



3 



space, subsequent recursion always stays within this subspace. This allows us significant 
reduction in rank of the required subspace and makes this method practically feasible. 
Further, this method with its terminator approximation, retains the Herglotz properties 
of the configuration averaged Green function. 



2. The augmented space recursion 



2.1. The recursion method 



The recursion method addresses inversions of infinite matrices jl]. Once a sparse 
representation of an operator in Hilbert space, if, is known in a countable basis, 
the recursion method obtains an alternative basis in which the operator becomes tri- 
diagonal. This basis and the representations of the operator in it are found recursively 
starting from = |_R) (8> |1) : 



where. 



\u2) = H\ui) - ai\ui) 

\Un+l) = H\Un) - an\Un) - (3l\Un-l) n>2 (l) 



= ~i — \ — ^ = 7 1 r (2) 



{Un\Un) {Un-l\Un-l) 

for n> 2. In this basis the representation of H is tridiagonal 



{Un\ 


\H\ 


\Un) 


{u 


n 


\Un) 



{Un-l\H\Un) / ifri \ n ^ , o 



and the averaged Green function can be written as a continued fraction : 



^GnRiE):^ = {ui\{Ei ~ H)-'\ui) 



E — ai 



PI 



E — a2 



PI 



E-aN-T{E) 

(3) 
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The asymptotic part of the continued fraction is obtained from the initial set of 
coefficients {a„,/5„} for n < N, using the idea of terminators For small values 

of we have large inaccuracies and the more the structure in the spectral density of H 
the larger is the N needed to maintain the window of accuracy required by us. 

2.2. The TB-LMTO Hamiltonian for random binary alloys 

In earlier communications we have described how to deal with random binary alloys 
P] within the framework of the tight-binding linearized muffin-tin orbital method (TB- 
LMTO) [Hlini and using the augmented space formalism (ASF). We refer the readers to 
the seminal papers mentioned earlier for details. Here we shall introduce only the salient 
features which will be required by us subsequently in this communication. We begin by 
setting up a muffin-tin potential with centres at the atomic sites i? on a lattice. Next, we 
inflate the muffin-tins into atomic spheres and start from a most-localized TB-LMTO 
representation of the Hamiltonian : 

H = ^ Crl Vrl + X! X! ^^rlSrl,r'L'^^r^l' Trl,r'L' 

RL RL R'L' 

The potential parameters, C and A, are diagonal matrices in the angular momentum 
indeces, and 

Crl =CtnR + CE (1 - ur) 

A](l= {Atf\R+{Aff\l-nR) (4) 

Here L is a composite label {i,m,ms} for the angular momentum quantum numbers. 
Ur is the random site-occupation variable which takes values or 1 with probability x 
or y {x + y = 1) respectively, depending upon whether the muffin-tin labelled by R is 
occupied by A or i?-type of atom. 

The Vrl and Trl,r'l' are projection and transfer operators in the Hilbert space H 
spanned by the tight-binding basis {\RL)}. 

2. 3. The augmented space formalism within the TB-LMTO 

We can associate the random variable hr with an operator Mr whose eigenvalues 
correspond to the observed values of ur, and whose corresponding eigenvectors, 
{|Oij), \1r)} are the state vectors of the variable. These state vectors of the set of 
random variables ur of rank 2^ span a space called configuration space $ with 
configurations of the kind | tt| ... |t ...) where 
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I Tr) = |0/?) + ^ \Ir) 

\iR) = ^/y |0r) - \Ir) 

If we define the configuration | || . . . | . . .) as the reference configuration, then any 
other configuration may be uniquely labelled by the cardinality sequence : {Rk}, which 
is the sequence of positions where we have a J, configuration. The cardinality sequence 
of the reference configuration is the null sequence {0} 

The augmented space theorem states that 

« AHur}) » = < {0}|A|{0} > (5) 



where, 

A = J... J A{{Xr}) l[dPi>^R) e $ 

P(Aij) is the spectral density of the self-adjoint operator Mr, which is such that the 
probability density of ur is given by : 

p{nR) = -- lim (tij | ({nR + i5)I - Mr) \ Tr) 

TT 0— ►0 ^ ' 

Mr is an operator in the space of configurations iJjr of the variable ur. This is of rank 
2 and is spanned by the states {| ti?), | Ir)}- 

Mr ^ X Vf + y vf + ^ (r^f + r^) (6) 

The expended Hamiltonian H in the augmented space is constructed by replacing all 
the random variables ur by the corresponding operators Mr. It is an operator in the 
augmented space ^ = 7i ® So the ASF maps a disordered Hamiltonian described 
in a Hilbert space 7i onto an ordered Hamiltonian in an enlarged space ^ , where 
the space ^ is constructed by augmenting the configuration space $ = Yl^ipR of the 
random variables of the disordered Hamiltonian together with the Hilbert space Ti of 
the disordered Hamiltonian. The configuration space $ is of rank 2^ if there are N 
muffin-tin spheres in the system. 

We may rewrite the expression for the configuration average of the Green operator as : 



<G(£;)> = <{0}|P + B + 'P|{0}> (7) 



where. 



A= Y.A[{E-Cl)//^l 



RL 



and 



Since, 



RL 
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c ^Y.'^[{e-Cl)/^l] VR®VL®{rfi+r^^} 
s 



X! X! Srl,r'L' Trri ® Tlu ® 1 

RL R'L> 



V^Y^^i^^T) Vr®Vl®T+Y.b[^T) Vr®Vl®VI 

RL RL 



RL 



A{V) = xVa + yVs i-e. the average of V 

B{V) = {1-2x){Va-Vb) 

C{V) = y^{VA - Vb) 



mn) - A{Ai'/')m)+c{Ai'/')\{R}) = |i} 



The ket |1} is not normahzed. We first write the above in terms of a normahzed ket |1) 
= [A(l/A)]"^/^ |1}. We now have : 



<G(E)> = {1\ (^E - A' - B' - C' - S') V) 
= (1| (e - H'f^y'\l) 



h'^^ = ^j^B' + C' + S' 



(8) 



where, 



A! = Y.{A{Cl/Al)/A{1/Al)} VR®VL®i 



RL 



-B' = E [i^L - E)/Al] I a (1/Ai)} Vr®Vl® Vl 



R 



RL 



c' = E [{.Cl - e)/Al] I a (i/A^)} Vr®Vl® {r« + T^f } 



RL 



= E E Srl,r'l' [A Trr, ® Tll' ® T 

RL R'L' 



(9) 
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This equation is now exactly in the form in which the recursion method may now 
be apphed. The computational burden is considerably reduced due to this diagonal 
formulation, the recursion now becomes energy dependent as is clear from the form of the 
effective Hamiltonian as shown in ((Tj) and discussed in pT| l22j. The recursion formalism 
of the ordered Hamiltonian was free of this constraint. This energy dependence 
makes the recursion technique computationally unsuitable because to obtain the Green 
functions we have to carry out recursion per energy point of interest. This problem 
has been tackled using seed recursion technique |221- The idea is to choose a few seed 
points across the energy spectrum uniformly, carry out recursion over those points and 
then interpolate the values of coefficients across the band. In this way one may reduce 
computation time. For example, if one is interested in an energy spectrum of 200 points, 
in the bare diagonal formulation recursion has to be carried out at all the 200 points but 
in the seed recursion technique one needs to perform recursions only at 10-15 points. 
The whole idea stems from the fact that in most of the cases of binary alloys, it is seen 
that the recursion coefficients and 6„ vary quite weakly across the energy spectrum. 
At this point we note that the above expression for the averaged -C Gll{E) 3> is exact. 

3. Symmetry reduction of the augmented space rank 

We mentioned earlier that recursion on the augmented space is not computationally 
feasible because of its large rank. For a binary alloy with N sites and with only s- 
orbitals on them, the rank of the augmented space is x 2^. Implementing recursion 
on this huge space for a sufficient number of steps to ensure accuracy is often not feasible 
on available computers. However, if we exploit the symmetry both of the underlying 
lattice in real space and of the configuration space (which arises due to the homogeneity 
of the disorder and arrangements of atoms on the underlying lattice), the rank of the 
irreducible part of the augmented space in which the recursion is effectively carried out 
becomes tractable. The conceptual advantages in ASF include apart from analyticity, 
translational and underlying lattice symmetries automatically built into the augmented 
space Hamiltonian. This allows us to involve the idea of utilizing symmetry operations 
present in both the real and configuration spaces, in the context of recursion method, 
reducing the rank of Hamiltonian drastically and making the implementation of ASF 
feasible. Since the augmented space recursion essentially retains all the properties of real 
space recursion but is described in a much enlarged space, it will be useful to consider 
symmetry operations in real space recursion first and then to consider those in the full 
augmented space. 

During the recursion, the basis member |m„) is generated from the starting state \ui) by 
repeated application of the Hamiltonian. If the starting state belongs to an irreducible 
sub-space of the Hilbert space then all subsequent states generated from recursion will 
belong to the same irreducible sub-space. Physically, we may understand this as follows: 
the recursion states carry information of distant environment of the starting state. 
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I 

Figure 1. Nearest neighbour shells around a central site on a square lattice. The 
colour coding indicates the symmetry weight of the sites : white represents sites with 
weight 1, gray with weight 4 and black with weights 8. 

If the Hamiltonian is nearest-neighbour only, then the state |m2), which arises by the 
apphcation of the Hamihonian on \ui) is a combination of states in the nearest shell with 
which \ui) couple via the Hamiltonian. Similarly, is a combination of n-th neighbour 
shell with which \ui) is coupled via the Hamiltonian. If 7^ is a point group symmetry of 
the Hamiltonian, then all n-th neighbour- shell states which are related to one another 
through the symmetry operator must have equal coupling to Hence, it is useful to 
consider among the n-th neighbour- shell states of which is a linear combination, 
only those belonging to the irreducible subspace and redefine the Hamiltonian operation. 

As an example, take a nearest-neighbour s-state Anderson model on a square lattice, 
with a binary distribution of its diagonal elements. We shall label the tight-binding basis 
by the position of the lattice points in units of the lattice constant : e.g. \{mn)) where 
m, n are integers 0, ±1, ±2, . . .. The starting state |(00)) belongs to the one dimensional 
representation of the point group of a square lattice. This state then couples with linear 
combinations of states on neighbour shells which are symmetric under square rotations : 

|(0m)} = (|(0m)) + |(mO)) + |(mO)) + |(0m))) /2, m>0 

|(mm)} = {\^{mm)) + \{fhm)) + |(mm)) + |(mm))^ /2, m > 

\{mn)} = (j{mn)) + \ {mn)) + \ {mn)) + \ {mn)) . . . 

+ \{nm)) + \ {nm)) + |(nm)) + |(nm))) /{2V2) 

m> 1, < n < m (10) 

The figure Q shows the grouping together of sites with the local symmetry of the square 
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lattice on the first three nearest neighbour shells of the central site. The first and second 
groups in (fTUj) are coloured gray in figure [T] and the last group coloured black. The labels 
on the groups (shown on the left sides of the equations (|Tn|l ) are confined only to the 
upper right quadrant of the lattice. 

If we go up to shells (for large N) there are about 2N'^ states in the diamond shaped 
nearest-neighbour cluster. However, there are only {N'^/A + N/2) ~ states with 

square symmetry. 

So within this reduction we can work only in 1/8-th of the lattice, provided we 
incorporate proper weights to the states to reproduce the correct matrix elements. If 
\R) and \R') are two states coupled to each other via the Hamiltonian, and both belong 
to the same irreducible subspace, and if |i?2), • • • , \Rwn) states obtained by 
operating on \R) by the symmetry group operations of the real space lattice, then Wr 
is called the weight associated with the state labelled by R. If we wish to retain only 
the states in the irreducible subspace and throw out the others and yet obtain the same 
results, we have to redefine the Hamiltonian matrix elements as follows : 



In the above prescription, the new irreducible basis introduced by us, reflects only the 
symmetry of the underlying lattice and holds good for a model system which has s-like 
orbitals only. But for a real system, the TB-LMTO minimal basis contains members 
with other symmetries as well. For example, in a cubic lattice with a spd minimal 
basis, we have basis members with s, p, Cg and t2g symmetries. The symmetry of 
the orbitals is reflected in the two-centered Slater-Koster integrals and this prohibits 
overlap integrals at certain positions, called symmetric positions with respect to the 
overlapping orbitals. A few of these symmetric positions for a simple cubic lattice are 
shown in figure El For the sake of clarity we have shown only the projections on x — y 
plane. We have indicated the positive and negative parts of the orbital lobes by different 
shades. From figure Efb) it is easy to argue that since the Hamiltonian is spherically 
symmetric and the product 0^(000)0^^(010) is positive in the upper right quadrant 
and negative in the upper left quadrant : {0, s\H\2,px) = 0. The same is true for 
{0,s\H\2,py) = 0. Similarly, {0,s\H\3,Pz) = {0,s\H\3,Pz) = 0. On the other hand, 
{0,s\H\l,p^) = -{0,s\H\l,p,) ^0. 

The above is an illustration, detailed arguments for the orbital symmetry has been 
discussed by Harrison fIU\ . 

The table ^ provides the conditions for obtaining the orbital based symmetry factor 
I3ll'{R — R') for a lattice with cubic symmetry. We shall describe the direction of 
R — R' by the direction cosines (imn). The table gives the details for the calculation of 
the symmetry factor PLLi{imn) on a cubic lattice. A look at the figures El and El shows 
the reason for introduction of this factor. 




(11) 
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Matrix 


Slater-Koster 


Symmetry 


Condition for it 


Element 


Expression 


Factor 


to be zero 


tss 


tssa 


r s^s 




ts,x 


f t 


Qlmn 

'~'^S,Px 


1= 


ts,xy 


^3 \m tsda 


olmn 
r's,xy 


1=0 or m=0 


^s^x^—y^ 




olmn 
1"^ s^x^—y^ 


l=m 




[n2-l/2(l2+m2)] tsd. 


Qlmn 
^'s,3z2-r2 


n= 1/V3 


^x,x 




fiimn 

r^X,X 


— 


^x,y 


hn (tppa tpp^^ 


olmn 
yx,y 


1=0 or m— 


^x^xy 


m [a/SI^ t„drr + fl-2P) t„dJ 


oimn 
1 x,xy 


m— 


t ^ ^1 ■y 

x,y^ 


\f% \mn (t„da -2 t^dw) 


olmn 
r-^x,yz 


1=0 or m=0 • • • 
• • • or n— 


^x^x^—y^ 


\ [v^/2 (P-m2) t^d. +(l-P+m2) t^d.\ 


olmn 
1 x,x'^—y'^ 


1=0 


^z,a;2— y2 




olmn 
r^z,x'^—y'^ 


n=0 ot \—m 


^a;,3z2— r2 


1 {[n2-(l2+m2)/2] t^da-^Tin^ W} 


olmn 
Px,3z^-r'2 


1=0 


trrii T.1I 


2>^m^ trfd<,+(l2+m2-4l2m2) tdd-K ■ ■ ■ 
■ ■ ■ tdd& 


olmn 
'Xy^xy 




^xy,yz 


\n{ ?)Vr? tdda^{^-^rn?\ tdd^} + (m'^-1) t,id/i 


olmn 
I xy^yz 


1=0 or n=0 


txy,x^—y'^ 


\m (l2-m2) { (3/2)idd^-2 i^^^ + (1/2) ^^^4 


d^"^\ 2 

"xy,x''—y'' 


1=0 or m=0 ■ ■ ■ 
■ ■ ■ or \=m 


tyz,x^—y^ 


mn{(3/2)(P-m2) tddA^+'^i^^-m^)] tddn ■■■ 
■ ■ ■+ [l+l/2(P-m2)] tdds} 


d^"^\ 2 

l^yz,x-'-y' 


m—0 or n—0 


txy,3z'^—r'^ 


V3 \m { [n2-l/2(l2+m2)] tdda-'^^ tddn ■ ■ ■ 
■ ■ ■ + l/2(lW) tdds] 


olmn 

r^xy,3z'^—r'^ 


1=0 or m—0 


^a;2— y2^a;2_j^2 


3/4(P-m2)2 t,,^ +[P+m2-(l2_^2)2] . . 
••• +[n2+l/4(P-m2)2] trfrf^ 


olmn 

t-'x^—y'^ ,a;2— 2/2 


- 


tx'^—y'^,3z^—r^ 


x/3 (P-m2){ [nV2-l/4(l2+m^)] t,,, -n^ t,,,- ■ ■ 
• • • +l/4(l+n2) tdds 


olmn 

/^a:2-j/2,322_r2 


l=m 


^32;2— r2,3z2— r-2 


[n'-l/2{Um')]' tdda +3n\l'+m') tdd.- ■ ■ 
■ ■ ■ +3/4(12+^2)2 


olmn 

/^322-r2,322_r2 





Table 1. Table showing the Slater-Koster parameters and the consequent symmetry 
factors The parameters for which the zero condition have been omitted can be 

obtained from the last column by permuting the direction cosines. 
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V 1 



Figure 2. sp-orbitals on a cubic lattice. Only projections on the xy plane are shown. 
The positive parts of the lobes are unshaded, while the negative parts of the lobes are 
shaded. The values of the symmetry factor for the overlaps {^Lvi^mn), defined in 
the text) are shown for each direction. The numbering of sites on the cubic lattice is 
shown in the inset of the right bottom. 

For the s-orbital overlaps the symmetry factor is always 1 (figure 1(a)). All six 
neighbours in the cubic lattice in the directions (100), (100), (010), (OlO), (001) and (OOl) 
are equivalent with the same overlap integrals tasa (only four are shown in the xy plane). 
We may then reduce the lattice and retain only the neighbour in the (100) direction and 
scale the overlap tssa by the appropriate weight (i.e. 2 in this case) as discussed earlier. 
For the sp^, however, the overlaps in the (010), (010), (001) and (001) are zero. This 
can be deduced immediately by noticing that the overlap products in these directions 
are positive for x > and negative for x < 0. On reduction, although they are related 
to one another by the symmetry operations of the cubic lattice, we cannot club the 
neighbours in the (100) and (100) direction with those in the (010), (010), (001) and (00 
1) and will have to put these overlaps to be zero by introducing the symmetry factor 
(5lu {imn) . 



{RL\H\R'L') ^ (3LL'{imn) {RL\H\R'L') (12) 

where j3LLi{imn) is if R' is a symmetric position of R with respect to L and L', 
otherwise it is 1. 
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Figure 3. Overlaps involving d orbitals in the eg and <2g symmetries. Only projections 
on the xy plane are shown. The positive parts of the lobes are unshaded, while the 
negative parts of the lobes are shaded. 



FigureElillustrates some of the overlaps involving the d states with Cg and t2g symmetries. 
With this simple reduction procedure the real space part can be reduced to 1/8-th, 
keeping only the sites in the (x > 0, ?/ > 0, 2; > 0) octant and suitably renormalizing the 
Hamiltonian matrix elements as described earlier. 

We still have not exhausted all the symmetries in the full augmented space. As 
discussed earlier, this space is a direct product of the real (lattice-orbital) space and 
the configuration space which are disjoint. As a consequence the symmetry operations 
apply independently to each of them. Since the disorder is homogeneous, the cardinality 
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Figure 4. Equivalent configurations on a cubic lattice. Configurations labelled by "f 
are shown as light spheres, while those labelled by [ are shown as dark spheres. The 
top twelve configurations are equivalent, while the bottom three are equivalent. The 
lower right inset shows the numbering scheme used for the lattice sites in the main 
text. 
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5 6 7 8 9 10 11 

Recursion Steps 

Figure 5. Convergence of the errors in energy moments of Ag^^^Pd^i^ as functions of 
the recursion step. 

sequence in configuration space itself lias the symmetry of the underlying lattice. To see 
this, let us look at the figure H] where we show a part of a cubic lattice where the central 
site is occupied by a configuration labelled by f, while two of the six nearest neighbours 
are occupied by configurations labelled by |, and four of them by js. We note that 
the twelve configurations in the first three rows of the figure, where the two |s sit at 
distances \/2 times the lattice constant, are related to one another by the symmetry 
operations of the cubic lattice. For example, the second to the fourth configurations 
on the top row can be obtained from the first by the rotations 7?.(7r/2, z), 7?.(7r, z) and 
7?,(37r/2,5) respectively. The configurations are described by cardinality sequences (as 
described earlier). The cardinality sequences for the four configurations on the top row 
of figure m are : {13}, {23}, {13} and {23}. From the figure it is easy to see that : 

{23} = 7^(7^/2, z) {13} ; {13} = 7^(7^, z) {13} and {23} = 7^(37^/2, z) {13} 

This equivalence of the configurations on the lattice is quite independent of the 
symmetries of the Hamiltonian in real space discussed earlier. Thus, in augmented 
space, equivalent states are \R®{C}) and the set |3fJi?®3?{C}) for all different symmetry 
operators 3? of the underlying lattice. Again the symmetry of the orbitals also rules out 
the operation of the Hamiltonian at certain symmetric positions discussed earlier. 

Once we have defined the reduced Hamiltonian, recursion on the reduced augmented 
space with starting state \ui) = \RL (g) {0}) gives the configuration averaged Green 
function directly. 

In order to give the readers a flavour of the reduction in storage and CPU time, we have 
carried out two sets of recursion calculations : First, a standard recursion on ordered 
Ag on a fee lattice, and second, an augmented space recursion with 15 seed points on a 
fee alloy Ag^oPd^o both with and without symmetry reduction. The calculations were 
done on a P4 machine with 1.13 GHz clock speed and 256 MB RAM. 
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Number of shells 



Figure 6. Showing how the number of lattice sites increase with increasing the number 
of shells in the real space map and its reduced version. This is for an ordered fee lattice. 



The first point to note is that convergence of the recursion technique is measured by 
the convergence of the energy moments of the corresponding density of states obtained 
from recursion 1^ • These energy moments are defined by : 



dE E'^niE) where / dE n{E) = 

-oo J —oo 

where Ug is the number of valence electrons. The convergence of these moments with 
the number of recursion steps N, from which the density of states n{E) and the Fermi 
energy Ep is calculated, is reflected in the errors 

AM„^ = |M„^-M„^-^| 

The figure El shows the convergence of the first three moments for a Ag^oPd^o disordered 
alloy on a fee lattice in a spd TB-LMTO minimal basis set. It is clear from the 
figure that for a convergence within a specified error window one has to carry out 
recursion for a specific large number of steps determined by the error window. Symmetry 
reduction allows us to carry out recursion eactly over much larger number of steps than is 
possible for ordinary recursion, given our computational resources. There was a similar 





Full lattice 


Reduced lattice 


System 


Shells 


Steps 


Sites 


CPU time 


Sites 


CPU time 


Model 


20 


20 


24739 


2.36 


3385 


0.54 


Ag 


20 


30 


24739 


47.15 


3385 


15.83 



Table 2. Comparison between system size and CPU time (in seconds) taken for 
recursion on a P4/256 machine for a full fee lattice and the reduced lattice in real 
space. 
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Number of shells 

Figure 7. Showing how the number of lattice sites increase with increasing the number 
of shells in an augmented space map on a disordered binary alloy on a fee lattice. 

observation made in an earlier estimate of the errors of the recursion method |24j. In 
order to carry out this many recursion steps we have to generate at least this many 
shells exactly around the starting site in order to avoid surface-like effects. 

For the case of the ordered fee lattice, the figure IHl shows how the size of the map increases 
as we increase the number of nearest neighbour shells from a starting site, both with 
and without reduction. Table 2 shows the details of the CPU time and storage space 
reduction for recursion after applying the symmetry reduction. We have carried out 
calculations both on a simple s-state system on a fee lattice , as well as for Ag (with 
spd minimal TB-LMTO basis) also on a fee lattice. 

For the calculation of a disordered binary alloy on a fee lattice, the figure [7| shows 
the enormous decrease in the size of the augmented space map after application of 
symmetry reduction for a seven nearest neighbour map on a fee lattice. TableEltabulates 
the reduction in storage and CPU time for a 7 shell, 11 step recursion in augmented 
space carried over 15 seed points using TB-LMTO potential parameters and structure 
matrix to built the Hamiltonian. The power of symmetry reduction of storage space 
is more evident in this example. Such reduction will allow us to stretch our nearest 
neighbour map up to 9-10 shells, stepping up our accuracy. The figure |H1 indicates the 





Full lattice 


Reduced lattice 


System 


Shells 


Steps 


Sites 


CPU time 


Sites 


CPU time 


Model (50-50) 


7 


7 


49476 


6.67 


6856 


1.21 


AgsoPdso 


7 


11 


49476 


1136.6 


6856 


336.27 



Table 3. Comparison between system size and time (in seconds) taken for recursion 
(using 15 seed points for AgsoPdso) on a P4/256 machine for a full fee lattice and the 
reduced lattice in augmented space. 
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Figure 8. CPU time taken per recursion step for a AgsoPdso disordered binary alloy 
on a fee lattice in a spd TB-LMTO minimal basis set. The times with and without 
symmetry reduction are shown. 

reduction in CPU time as we increase the number of recursion steps. Further, since 
the number of sites in the map decrease, the number of individual apphcations of the 
9x9 Hamiltonian also decrease significantly, as do the number of operations involved 
in taking various inner products. This will reduce the inherent cumulative error of the 
recursion technique and lessen the probability of the appearance of ghost bands which 
often plague recursion calculations. 

4. Conclusions 

We have shown that recursion calculations can be carried out much faster and for many 
more recursive steps exactly, if we perform the recursion on a subspace of the original 
augmented space reduced by using the symmetries of both the underlying lattice and 
random configurations on the lattice. This will allow us to obtain results for disordered 
binary alloys with enhanced accuracy required for first-principles, self-consistent, density 
functional based calculations. In this communication we have described the details of 
the implementation of this symmetry reduction and the modifications required in the 
standard recursion method. We propose using the symmetry reduced version of the 
augmented space recursion in our future work on disordered alloys. 

References 

[1] Saha T., Dasgupta I and Mookerjee A., 1994, J. Phys.: Condens. Matter 6 L245 
[2] Saha T., Dasgupta I. and Mookerjee A., 1996 J. Phys.: Condens. Matter 8 1979 
[3] Mookerjee A., 1973 J. Phys. C: Solid State Phys. 6 L205 ; 1973 J. Phys. C: Solid State Phys. 6 
1340 

[4] Haydock R., Heine V. and Kelly M.J., 1972 J. Phys. C: Solid State Phys. 5 2845 




Symmetry reduction in the TB-LMTO-ASR for random binary alloys 



18 



[5] Andersen O. K., 1975 Phys. Rev. B 12 3060 

[6] Jepsen O. and Andersen O. K., 1971 Solid State Commun. 9 1763 

[7] Mookerjee A. and Prasad R., 1993 Phys. Rev. B 48 17724 

[8] Saha T. , Dasgupta I. and Mookerjee A., 1994 Phys. Rev. B 50 13267 

[9] Saha T. and Mookerjee A., 1995 J. Phys.: Condens. Matter 8 2915 

[10] Alam A. and Mookerjee A., 2003 Phys. Rev. B (to be published) ; cond-mat/0307697" 

[11] Beer N. and Pettifor D. G. (1982), in The Electronic Structure of Complex Systems, ed. P. 

Phariseau and W. M. Temmerman, NATO ASI Series B, v.ll3, p 769 

[12] Lucini M.U. and Nex C.M.M., 1987, J. Phys. C: Solid State Phys. 20 3125 

[13] Kaplan T. and Gray L.J., 1976 J. Phys. C: Solid State Phys. 9 L303, L483 

[14] Kaplan T. and Gray L.J., 1977 Phys. Rev. B 15 3260 

[15] Ghosh S., Leath P.L., and Cohen M. H., 2002 Phys. Rev. B 66, 214206 

[16] Kumar V., Mookerjee A. and Srivastava V.K., 1982, J. Phys. C: Solid State Phys. 15 1939 

[17] Srivastava V.K., Choudhry V. and Mookerjee A., 1983 J. Phys. C: Solid State Phys. 16 4555 

[18] Haydock R. and Te R.L., 1995 Computer Phys. Communication 90 81 

[19] Gallagher J. (1978), Ph. D. Thesis, University of Cambridge, U. K. 

[20] Harrison W.A., 1999 Elementary Electronic Structure (World Scientific, Singapore) p 546. 

[21] Biswas P. Sanyal B., Mookerjee A., Huda A., Haider A. and Ahmed M., 1997 Int. J. Mod. Phys. 
B 11 3715 

[22] Ghosh S., Das N. and Mookerjee A., 1999 Int. J. Mod. Phys. B 21 723 

[23] Haydock R., Solid State Physics ed. H. Ehrenreich, F. Seitz and D. TurnbuU (Academic Press, 
New York) (1980) 

[24] Chakrabarti A. and Mookerjee A., 2002 J. Phys.: Condens. Matter 14 3211 



